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Abstract 

Weyl functions conveniently describe the evolution of wave coherences in periodic or quadratic 
potentials. In this work we use Weyl functions to study the "Talbot-Lau effect" in a time-domain 
matter-wave interferometer. A "displacement diagram" is introduced to analyze and calculate the 
matter-wave interference for an atomic cloud in a quadratic potential that interacts with a sequence 
of short optical standing wave pulses producing an atomic grating echo. Unlike previous treatments, 
this new approach allows the atomic ensemble to have an arbitrary initial phase-space distribution, 
and the standing wave grating vectors to span three dimensions. Several examples are discussed to 
illustrate the convenience of the diagrammatic technique including the following: a two-dimensional 
Talbot-Lau effect, the shift in the echo time and the recoil phase for the interferometer perturbed 
by a quadratic potential; and the realization of a time-domain "Lau effect" using a pulsed harmonic 
potential. The diagrammatic technique is applicable to diffraction gratings with arbitrary grating 
transmission functions. We conclude the paper with a general discussion on the Weyl function 
representations of matter-wave coherence, and relate the conservation of matter-wave coherence 
with the conservation of purity ? = Tr{p'^) that distinguishes decoherence effects from dephasing 
effects. 

PACS numbers: 39.20+q, 03.75.dg 
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I. INTRODUCTION 

In quantum phase space methods, operators in a Hilbert space are mapped to functions 
in a "mock" phase space [l| . One example is the Wigner representation which is particularly 
useful for establishing quantum-classical correspondence. In the Wigner representation, the 
Wigner function represents a density matrix operator p with a quasi-probability distribution 
in classical phase space. The Wigner functions are conveniently applied in semiclassical 
methods j^, sl, to visualize or characterize quantum states, and experimentally allows one 
to recover information of quantum systems through tomographic measurements of their 
phase-space distributions jj]. 

In the related Weyl function representation [l, 5], the Hilbert space operators are mapped 
to functions in a "reciprocal phase space" where each function is simply a Fourier transform 



of those in the Wigner representation. In particular, Weyl function iy(q, x) is t 
transform of the Wigner function w(r, p) over the phase space coordinates (r, p) 



le Fourier 
The 



g. 
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Weyl function can also be expressed as: 

PF(q,x) = Tr(pD(q,x)), 
£>(q,x) = e^('i-'-^P). 

Here r and p are the position and momentum operators, and D(q, x) is the displacement 
operator js], that shifts the density matrix by /iq in momentum and x in position. Thus the 
Weyl function iy(q, x) is a phase-space correlation function that measures the coherence 
of matter-waves. Though both Weyl and Wigner functions completely specify any density 
matrix operators, application of Weyl function, or the Fourier transform of the Wigner 
function, has been much less discussed before 

In this work, we demonstrate the power of the Weyl function technique by applying it 
to an important class of wave interference phenomena, the Talbot-Lau effects associated 
with matter-waves in quadratic or periodic potentials. Various exarnples of Talbot-Lau type 
interference, including those closely related to recent experiments 8|, l9|, llOl] , are discussed. 
To our knowledge, this work gives the first vectorial and diagrammatic formula of Talbot- 
Lau effects. Apart from the expected application in matter-wave optics, the formula may 



also provide a convenient framework for the design of various lens-less imagine systems 



A. Talbot-Lau effects 

Periodic light wavefronts can repeatedly reconstruct themselves during free propagation. 
The phenomena is known as the Talbot effect that was discovered by Henry F. Talbot [l^ in 
1836. A demonstration of the Talbot effect usually starts with a collimated monochromatic 
light with wavelength A. The light passes through a transmission grating with a grating con- 
stant d to generate a periodic wavefront. Full or partial re-constructions of the transmission 
wavefront at integer or fractional multiples of the Talbot distance = (P/X are observed 
downstream. Many intricate and elegant features of the Talbot phenomena have been dis- 
cussed in the literature Instead of using a collimated or spatially coherent input, in 
1948 Ernst Lau [l4] found that spatially incoherent light, propagating through two gratings 
with same grating constant d separated by a distance Zi2, can generate an interference pat- 
tern at the focal plane of a lens (Fraunhofer diffraction) placed downstream from the second 
grating. The contrast of the interference is maximized if takes integer multiples of zt/2. 
More general than Fraunhofer diffraction, in the Lau effect the interference pattern appears 



at a distance 2:23 = p/qzi2 downstream from the second grating where p, q are two integers. 
Lau's original discovery corresponds to p/q = N ^ 00 14. |]^. 
Instead of using hght waves, both the Talbot Q, ll 



llTj l and the Lau effects 18| have been 
demonstrated using matter-waves and in particular using neutral atoms. Interferometry 
with multiple gratings using the Lau effect is referred to as Talbot-Lau or Talbot-vonLau 
Interferometry (TLI). The Lau effect is of particular interest for observing matter- wave 
interference in the circumstances where a bright and coUimated matter-wave input is not 



available. Ta^ 
molecules 16 



bot-Lau inter: 
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of large molecules 



20 



Urometers are used to demonstrate wave properties of atoms and 



211 ]. to study decoherence [22, |23| and to measure the polarizability 



2J] . Talbot-Lau interferometers have also been explored to sense inertial 



forces with matter-waves 



25 



Based on Fresnel diffraction, matter-wave Talbot and Lau effects can also happen in the 
time domain instead of in the spatial domain, due to the equivalence of Schodinger's equation 
and the paraxial wave equations. Matter-wave Talbot 17|| and Lau 261] effects in the time 
domain have been observed with atoms cooled down to a thermal velocity of centimeters per 
second or lower. In the time domain, spatially separated gratings are replaced by successive 
pulses of optical standing waves acting as either phase or amplitude gratings for the matter- 
waves. Time-domain interferometry using optical standing wave pulses with laser-cooled 



atoms were pioneered by the authors of ref. 



261 ]. In that experiment, off- resonant standing 



wave pulses act as thin phase gratings to an atomic cloud. After each standing wave pulse 
the atoms bunch to the potential minima of the standing wave to form a density grating. 
The density grating washes out due to velocity spreading of atoms. However, with two 
successive standing wave pulses separated by time T12 applied to atoms, the atomic density 
grating revives around time t = T12 + T23 if T23 = T12 and more generally if T23 = NTu, 
where N is an integer. The contrast of the revived density grating is a periodic function of 
T12, similar to the spatial Lau effect. The revived atomic density grating is probed using a 
Bragg-scattering of a probe light to induce a grating echo 28|, and the interferometry based 
on this setup is referred as a time-domain grating echo interferometer {2^. Experiments 



and the theory of the interferometry were further developed in ref. [29|] that generalizes the 
technique to include multiple standing wave pulses. Recently a direct spatial image of sub- 
wavelength spatial interference pattern due to the Lau effect after two standing wave pulses 
(act as either amplitude or phase gratings) is achieved using an on-resonant standing wave 
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mask [30|, l31 |. 



Developed theories of matter-wave Talbot-Lau effect (TLE) usually simplify the analy- 
sis by making particular assumptions on the coherence properties of input matter-waves. 



Closely related to the spatia! 
ings, in the first approach 32 



(beam) Talbot-Lau setup with amplitude transmission grat- 



33 



34l | the atom flux through each opening of the first grating 
are treated essentially as individual point sources, from which the Fresnel diffraction on the 
image plane is calculated. An incoherent sum of the contribution from all the openings 
at the 1*** grating gives the final interference pattern. Instead of incoherent summation of 



point sources, the second approach [19|, l29|, |35| assumes the input atoms as an incoherent 



sum of plane waves. However, a general theory of Talbot-Lau interferometry applicable to 
matter-waves with arbitrary initial spatial coherence hasn't been presented in the literature. 

The approach to TLE in this work is motivated by an obvious link between the Weyl 
function [Eq. ([1])] and the time-domain matter- wave TLE due to a sequence of standing wave 
pulses. In the following we consider a single atom picture with atomic motion characterized 
by a density matrix operator p. TLE corresponds to a revival of atomic density grating at 
a certain time t. Since the revived density grating is spatially periodic, we shall consider its 
particular Q component given by: 

PQ = Tr(p(t)e'Q-^) = W^(Q,0,t), (2) 

where iy(q, x, t) is the Weyl function defined in Eq. ([1]). 

To express Eq. ([2]) in terms of an arbitrary initial condition described by p(0), one can 
propagate /3(0) forward in time by considering the interaction of atoms with a sequence of 
standing wave pulses. Alternatively, we shall define the time-dependent displacement op- 
erator Z)(q, X, t) = e**^'' ''*^*)~^ P^*^^ in the Heisenberg picture and propagate the displacement 
operator D(Q,0,t) backward in time. Here we take the latter approach which is advan- 
tageous for two reasons: 1) It applies to general initial conditions. 2) As we will see, the 
equation of motion for -D(q, x, t) due to pulsed periodic potentials, during free-propagations 
and due to a quadratic potential can be integrated easily. Since our analysis is generally 
applied to different types of interferometers using TLE, we shall generally consider an inter- 
ferometry configuration using TLE as a Talbot-Lau Interferometer, that will be referred to 
with the abbreviation TLI in this work. 
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FIG. 1: Schematic of an N-pulse grating echo Talbot-Lau interferometer 
B. Outline 

In what follows the paper is organized into three sections. In the first section we introduce 
a diagrammatic technique to calculate the grating echo of an N-pulse TLI in free-space. In 
the second part we discuss the application of the theory to matter-waves in a quadratic 
potential. In the discussion section we extend the formula to include arbitrary grating 
transmission functions, discuss the Weyl function representation of matter-wave coherence, 



and compare the approach here with other diagrammatic techniques 36 1. 



II. GRATING ECHO DUE TO TLE IN FREE SPACE 
A. The model 

In this section we consider TLE due to the simplest possible interaction between 2- 
level atoms and a sequence of off-resonant standing wave light pulses. The scenario is 
summarized in Fig. [TJ an atomic sample localized at the waists of multiple standing wave 
fields is subjected to a sequence of N-1 standing wave pulses. At time T/v a grating echo is 
induced with the amplitude proportional to a particular spatial Fourier component of the 
revived atomic density grating. The assumptions in the model are summarized below. For 
notational convenience, we set ^ = m = 1. 

Al. 2-level atoms are always in their adiabatic ground state, with any spontaneous light 



scattering ignored. Standing waves pulse i with the time-dependent light shift 37 1 



potential = Qi(t) cos(Qi-r + y9j) acts as a phase grating to atoms. Here Qj = k" — k| 



b 



is the k-vector of the standing wave. k°, are the k-vectors of the travehng wave hght 
fields, whose electric field amplitudes are referred to as i?? and E\. And (fi determines 
the position of the standing wave nodes in the lab frame. We assume the light fields 
have Gaussian beam profiles and the atomic sample is located at the center of each 
beam. The cross-sections of all the beams are characterized by a waist size Wq, that is 
much larger than the size of the atomic sample. We also assume all the and k^ to 
have significant differences in propagation directions such that any linear combinations 
of QjS have magnitude either zero or on the order of with A the optical wavelength. 

A2. Standing wave pulses are in the Raman-Nath regime, e.g., the positions of atoms are 
considered frozen during each pulse duration. We thus can use Dirac delta functions 
to effectively describe the temporal profile of the pulses, e.g., VLiit) = 9i6(t — Tj) that 
is characterized by its strength 9i and the time of arrival Tj. 9i is referred as the pulse 
area. 

A3. The atomic ensemble is assumed to be dilute such that atom-atom interactions are 
ignored. The ensemble is assumed to be dilute in phase space as well so that the single 
atom picture is applied without conflicting with quantum statistics. We can thus use a 
Dirac ket {ip) to label the matter-wave state, that obeys the single atom Schrodinger's 
equation. 



Given a fixed amount of atoms participating in the interaction, the spatial density 
distribution of the atomic ensemble is given by p(r,t) = NatomPr,r(t), where pr,r'if) = 
('?/'* (r, t)ip{r', t)) is the single atom density matrix for the atomic ensemble in real space. 
Impacts of !fi on the matter-wave dynamics require only relatively simple algebra, so 
we may drop ipi during the derivation and recover their values if necessary. 

A4. A grating echo refers to the Bragg scattering of probe light pulse (the last graph in 
Fig. [1]) off the revived atomic density grating. The amplitude of the grating echo 
is evaluated under first order Born approximation. We assume the evolution of the 
atomic density grating is not perturbed by the scattered light during the probing 
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interval. Using a slowly varying amplitude approximation, the coupling between Ej^ 
and the probe light field E% around time T/v is governed by the equation: 

{dt - cd.jE'}, = iauNatomPrAty'^''"'Etf. (4) 

Here xq measures the distance along the direction Qn, and a is the dipole polariz- 
ability of ground state atoms at standing wave frequency u. If we ignore the light 
propagation retardation, the overlap between the Bragg-scattered field with the E'^ 
can be written as: 

E'ATn) = -^V^^^E%iT^) J dhp{r, T^)e^'^--\ (5) 

where rj is a. geometric factor of order unity depending on the overlap between the 
Gaussian mode Efj and E'}^. E'^^iT^) in Eq. ([5]) gives the magnitude of the grating 
echo, and is considered to be the output of the TLI. 

B. The free space displacement diagram 

1. The displacement operator and the Weyl function 

We define g = rj^^^^^^^Ef^iT]^). Recalling Eq. ([T]), we define the displacement operator 
L'(q,x) = e'^'^-^-''-^^ and let iy(q,x,t) = Tr(p(t)l)(q, x)), so that Eq. ^ can be rewritten 
in a compact form using the Weyl function: 

E'k{TN) = -tgW{ClN,Q,TN). (6) 

To calculate the TLI output Elj{T]\f) , we might integrate the matter- wave equation 
Eq. ([3]) to express p{T j^) in terms of p(T^). Here instead we shall take the Heisenberg picture 
and integrate the equation of motion for the displacement operator £)(Q,x) (Appendix A) 
so that iy(QAr, 0, T/v) can be expressed in terms of iy(q, x, Ti). During free-propagations 
we have, 

W(q,x,t) = l^(q,x-qt,0). (7) 

With the notation 7^+^" = ± 0+ to remove ambiguities due to delta functions, trans- 
formation of iy(q, X, t) due to a standing wave pulse is written as: 
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FIG. 2: Basic displacement diagrams: a) A propagating displacement line [Eq. ([7])]. b) A scattering 
vertex with an incoming and outgoing displacement line, the scattering amplitude is evaluated at 
the vertex position [Eq. ([8])] . The dashed arrow represents the action of a standing wave diffraction, 
c) A probing vertex that stops an input displacement line [Eq. The echo amplitude can be 

evaluated at the crossing between the incoming displacement line and the x-axis. The solid arrow 
represents the Bragg-scattering (grating echo) process. The direction of the dashed line array 
indicates the slope of the displacement line x(t) = Qt. 



(q, X, T+) = J2 sin ^^)W^(q - nQ„ x, Tr). (8) 

2. Basic displacement diagrams 

We shall consider Eqs. ([7]) and ([HD in a space (q, x) that may be called the reciprocal phase 
space. From Eq. ([7j), free propagation of W{q^,x.) is simply a shearing in reciprocal phase 
space along x. From Eq. ([8]), a standing wave pulse "scatters" iy(q, x, t) into multiple 
duplicates shifted along the q coordinate by nQ weighted by J„(2^sin^). We define 
W {qn{t) , yin(t)) = {W{q^,yi,t))n, e.g., the component that is shifted by nQ due to standing 
wave scatterings, and represent its dynamics using a point (q(t),x(t)) that is moving and 
scattered in (q, x) space. The trajectory can be projected onto x coordinates in a x — t 
diagram, where the slope of x(t) represents the magnitude of q(t). Equation ([7]) and ([8]) are 
diagrammatically summarized in Fig. [2](a) and (b) as the propagation diagram {displacement 
line) and scattering diagram {scattering vertex) for the Weyl Functions. In addition, to 
calculate i?^(T/v) fully diagrammatically, we introduce the probing diagram {probing vertex) 
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FIG. 3: A diagram corresponding to a particular displacement line in a general N-pulse TLI. 
in Fig. [2](c) to represent the Bragg-scattering given by Eq. (Q. 

3. The displacement diagram for an N-pulse TLI 

With further explanations coming shortly, we sketch a representative path of (q(t),x(t)) 
that contributes to the grating echo in the setup sketched in Fig. [TJ The path starts from 
(qi,xi) at time t = Tf, propagate toward (q^r = Qjv,xjv = 0) at T/v, scattered by N-1 
scattering vertex at Xj = x(rj) with pulse area 6i and standing wave k- vector Qj, with 
diffraction orders ni,i = 1, N — 1 assigned at each scattering vertex, with njv = 1 and 
Tij = Ti — Tj. By evaluating the Bessel functions at each vertex using the rules in Fig. [2l 
the contribution of the path to TLI output is given by: 

E%{Tj,) = -ig E [n J„.(2ft,sinto)]H/(q^,Xi,t = Tf), 

{n^} t=l (a\ 
N N V') 

with : qi = E %Qj,Xj = E njQjTj^i. 

j=i j=i 

In Eq. ([9]) we include a summation of all possible {nj} paths that connect (qi, xi) at t = T{~ 
with (qAT = Qiy,xiy = 0) at T/v. 

Equation Qj can be evaluated with the displacement diagram in Fig. [3] and the rules 
given in Fig. [2l This is the first major result of this work. To proceed to a calculation 
based on any initial atomic states, we still need to specify a particular (qi,Xi) at time T^^. 
Clearly, iy(qi,xi) needs to be non-zero for a non-vanishing contribution to Eq. ([9]). In the 
next subsection, W{qi,:Ki,t = T{~) for a "uniform and broad" atomic sample is studied and 
we find that (qi = 0, Xi ^ 0) is required. 
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FIG. 4: Diagram restricting the displacement line at the input of the first scattering vertex [Eq. (jl2l) 
and Eq. (jl4p ]. The blue Gaussian on the left represents the matter- wave spatial correlation function 
due to Maxwell's velocity distribution with coherent length Ac . Here qi = and xi need to be 
within the coherent length of the atomic sample. 

C. Momentum conservation and echo condition 

1. Input displacement line 

This subsection considers the possible values for (qi,xi) that determine the input dis- 
placement line in Fig. [3l We shall consider a "uniform and broad" atomic sample, specified 
by the phase-space distribution with the Wigner function: 



where g{r) is a spatial distribution function, assumed to be "uniform and broad", e.g., 
uniformly distributed over many optical wavelengths A. This requirement is more general 
than simply requiring a distribution over a large spatial extension, /(p) gives a velocity 
distribution. 

The Weyl function evaluated at (qi,xi) is given by: 



Here qi can only have discrete values specified by Eq. ([9]), which is either zero or much 
larger than the spatial frequency of g{v) according to the assumption (Al) in section Hi A[ 
A Non-vanishing integration in Eq. (fTTl) over r thus requires the condition of momentum 



Wi3{y,y>) = ^(r-/?p)/(p). 



(10) 



iy(qi,xi,t = Tf ) =< e'li-^-^i P > 
= / (i3rrf3p«;^(r, p)e^('ii '-^i P). 




(11) 
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FIG. 5: Echo line in Fig. [3l Here the black, blue and red vertical bars indicate the arriving time 
of the first pulse Ti, the echo time T'^'^^°, and the arrival time of the A^*'^ pulse T/v respectively. 

conservation: 

N 

qi = E^.Q. = o- (12) 

The restriction from Eq. fll2l) simplifies Eq. fill I) and we end up with: 

W^(qi,Xi,t = rr) = /rf3p/(p)e-iP 

= F(xi) (13) 



where F{x.) is the Fourier transform of /(p). The velocity distribution function /(p) specifies 
a coherence length Ac. As a simple example, in the third line of Eq. f|T3l) we consider an 
isotropic Maxwell distribution of atom velocity as /(v) oc e 2^ , and we define the coherence 

N 

length Ac = — . Non- vanishing value of Eq. f|T3l) requires xi = ^jQjTj,i to be on the 
order of Ac. Combining with Eq. (fT^ . we have: 



N 

xi = ^n,Q,T,,i~o(Ac). (14) 

J=2 



Equations f[T^ and f[T^ require the input displacement line in Fig. [3] to start with (qi = 
0,xi ~ o(Ac)) in reciprocal phase space such that iy(q, x, Tf) is non-zero and can be 
scattered toward (QAr,0) at time Tjy. This is summarized in Fig. |4] as a supplementary 
diagram to Fig. [31 With the basic diagrams in Fig. [2] and with the input displacement 
line specified in Fig. HI output of an arbitrary N-pulse TLI can be evaluated using the 
displacement diagram in Fig. [3l 
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The condition of xi ~ o(Ac) in Eq. (HM is different for specific initial condition with 
different temperatures, or initial atom velocity spreadings. In the next subsection, this 
initial-condition-dependent requirement will be concretized with an "echo condition" for the 
displacement diagram. 

2. Echo condition and echo time 

In the displacement diagram given by Fig. [3l the displacement line connects the point 
(qi = 0, xi ~ 0) at t = Tf with the point {q^ = Qn, xat = 0) at T/v. If Tj , z = 1, — 1 
are fixed, then xi = xi(TAr) is a function of the detecting-pulse arriving time T^. 

With a particular TLI setting described by the parameters {(njQj,Tj),z = 1,...,A^}, if 
xi defined in Eq. (Q is along the direction of Qat, then by varying the time T/v, xi can be 
arbitrarily close to zero. The echo condition is defined as follows: In a general TLI setup 
[using a "uniform and broad" atomic sample specified by Eq. ([T0|) ] with the first N-1 pulse 
parameters set to be {njQj, Ti,i = 1,...,A^ — 1}, a revival of atomic density grating with the 
k-vector Qat requires the existence of an echo time T^'^'*° such that 

N 

xi(Tiv = = njQjTj^i = 0. (15) 

i=2 

In the displacement diagram, an echo line is a specific displacement line that connects 
(qi = 0,xi = 0) and (qiv = Qn,^n = 0). For the echo line, the input Weyl function 
evaluated by Eq. f|T3l) is maximized. An echo line for the N-pulse diagram in Fig. [3] is given 
by the dotted line in Fig [51 

3. Envelop function of grating echo 

With the echo time defined by Eq. (fT5l) and with the momentum conservation condition 
Eq. ( |T2l) . it is usually more convenient to write: 

xi(r) = Qnt, (16) 

with r = Tjv — T'^'^^° = We can then express the interferometer output as a function 

of r , the time offset of the probing time T/v from the echo time T^'^^°. From Eq. we 
see the grating echo given by Eq. ([9]) has a r dependent envelop function F{Q]^t) with its 
magnitude a symmetric function with respect to r = and reaches unity at r = 0. 
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4- Total interrogation time 



Notice that when the echo condition Eq. (fTS!) is satisfied for {Ti = 0, T2, ...,Tn = T^'^^° = 
T}, it must also be satisfied for T'{si = 0, S2, s^r = s'^'^^° = 1}; with Sj oc Tj for arbitrary 
T'. T shall be called the total interrogation time of the TLI. 

5. Vanishing grating echo at the echo time 

In Eq. ([9]) the only initial condition dependent term is the last term in Eq. iQ. For 
T/v = T'^'^^° such that W{qi = 0,xi = 0) = 1, the TLI output is not only maximized, but 
also becomes initial-condition independent. This is the feature of a "white light interference" 
as the interference is independent of the degree of input coherence. However, as noticed by 
the authors of Ref. 29|, if is set to be exact the echo time T'^'^^°, the TLI output given 
by Eq. ([9]) vanishes since in the pre-factors of Eq. ([9]) J„^(2^i sin ^^^) becomes zero for 
non-zero rii. Thus TLI using phase gratings are not "perfectly- white-light" interferometers. 

The vanishing grating echo at the echo time is understandable, as we notice the revived 
fringe at time T/v is a lens-less image of atomic density grating at time t = T/v — T^^'^° by 
the following N-2 gratings. If T/v is equal to T'^'^^°, the atomic distribution to be imaged is 
at t = Ti, the first pulse arriving time. Since the first grating is a phase grating, the input 
matter-waves only acquire a periodic phase at t = Ti' but keeps a uniform density profile. A 
density grating is developed only after a 6'1-dependent time during which the atoms bunch 
to the nodes of standing wave. The vanishing of the density grating at the echo time is a 
feature of TLE using phase gratings. Extension of interferometry with phase gratings to 
those with amplitude gratings will be discussed in the last section of this paper. 

For a hot atomic sample with Ac << A, we may Taylor expand the J„^ term in Eq. (!20|) 
to have cx t'^^F{Qi^t). For rii = 1 and a Maxwell velocity distribution of /(p), we get 
a dispersive-s hap ed backscattering curve as a function of r, consistent with experimental 



observations 26 1. 



6. Minimum number of pulses in D dimensions 

Finally, we notice that Eqs. (fT2|) and (fT5|) are two vector equations. In a particular 
experimental setup with fixed grating vectors Qj and timing Tj for each standing wave pulses, 
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a non-vanishing grating echo requires Eq. (fT2|) and Eq. (fT5|) have non-zero integer solutions 
for {ni,i = 1,...,N — 1}. Generally, for {Q^,^ = 1,...,A^} which span a D dimension, a 
minimum number of standing wave pulses are required to satisfy both equations, given by 

N^UD) =D + 2. (17) 

For > A^rnin; uiore than one set of integer solutions for rii may satisfy both equations. The 
contributions from different pathways coherently add up and may interfere. An example 
illustrating this interference effect (a "balanced 4-pulse" configuration) will be given later 
in this section. 



D. Interferometry observable 

1. Recoil phase 

N-l 

We now come to the pre- factor " n Jn, (2^i sin ^)" in Eq. (Q. We shall define: 

i=l 

irecoil Qi ' 



(18) 



as the recoil phase at each scattering vertex in Fig. [31 referring to the cumulative phase shift 
due to the kinetic energy of recoil effects. Clearly, 0™^°*' follows Xi and increases linearly 
with T, so its slope gives a temporal periodicity of the grating echo amplitude in Eq. ([9]). 
The periodicity of the gratin g ec ho amplitude was measured experimentally to determine 



the recoil frequency of atoms 
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38 



39] 



2. Acceleration induced phase 

We can generalize Eq. ([9]) in a constantly accelerating frame (Appendix A). In this case 
Eq. ([7]) is modified to include a phase factor: 

(q, X, t) = W{q, X - qt, 0)e'f^''-'^>^\ (19) 

where a is the acceleration constant. Equation is modified accordingly; we have, 

E'^{T^) = -zgY.[I[ J„.(20.sin^^)]iy(qi,Xi,Tf)e ' » , (20) 
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where x(t) is given by the displacement hne trajectory in the x — t diagram as shown in 
Fig. [31 In Eq. fl20l) we also include the phases ipi for each standing wave pulse explicitly. 
Equation fl20p indicates that the displacement line x(t) in Fig. [3] can be interpreted as the 
relative displacement between two diffraction path pair - external perturbations shift the 
relative phase between the two paths the same way as a Mach-Zehnder interferometer. 

3. Rotation induced phase or Sagnac phase 

Due to the equivalent principle, it is always possible to choose an inertial frame to describe 
the atomic motion, and attribute the inertial forces to the motion of the standing wave field 
classically, as long as the space-time curvature is negligible over the extension of the atomic 
ensemble. With this in mind, consider Fig. [1] with grating vectors rotated at an angular 
velocity fl along a rotation axis across the atomic sample. In this frame, the standing wave 
k-vectors Qi(Ti) is rotated relative to QiiTi) by an amount (5Qj = f2 x QjTi^i. Accordingly, 
Eq. (fT2l) is modified to be 

(21) 

5x1 (T^ = T-'^") = E X Q,T2 . 

In Eq. fl2Tl) the first line is zero due to Eq. f|T5|) . We end up with a non-zero displacement 
xi, which can be put together into Eq. (ITT!) while assuming a mean velocity of the atomic 
ensemble to be < v >. We end up with the rotation induced phase shift for the phase factor 
in Eq. (|20|): 

N 

= (< V > xfl) • ^ ^jQ^TI^. (22) 

i=i 

Equation ( l22l) may also be derived in a rotating frame and we consider the atomic Coriolis 
forces. On the other hand, equation f[T^ can be derived the same way as Eq. fl^2]) by 
considering the variation of the standing wave phases due to accelerations. 

E. Examples with a single grating vector 

In this subsection we give three examples of displacement diagrams for calculating TLI 
outputs. The most common experimental setup is with N standing wave pulses from a single 
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TLI name 



equal time 
3-pulse 



balanced 
4-pulse (n) 



elongated 
4-pulse 



Ref. Fig. Fig. El Fig. Fig. [8] 



J_l(29i sincjQT) J2re-3(26l sinUQT) J_ i (261 sin uqt) 

^ ' 72(292 sin^QCf +t)) j^_3„ (262 sin(a.Q { i^Zi^^Z + t))) Ji (292 sin a;Q (T, + t)) 

J„{2e3sin(ijQ{5 Ji(2e3sina.Q{T, +t)) 



97Q^? 6^r2 Ts{T-Ts) 

TABLE I: Interferometry setup name (first row), reference figure (second row), amplitude (third 
row) and acceleration induced phase shifts (fourth row) of the pq revivals in Figs.[6l[71[8l vq = 

~ '^T' '^c = ^ and u is the average thermal velocity of an atomic sample. The amplitudes 
of the revival is a product of all the Bessel functions in each column. In the third column the 
"balanced 4-pulse" has multiple loop contributions labeled with index n - the total amplitude is a 
coherent summation over n contributions including the phases, q (fifth row) gives the modification 
of revived k- vector due to a weak quadratic perturbation V = ^moofx'^. 

standing wave field. In this case there is only a single grating vector, e.g. Qj = Qbx. From 
Eq. (1171) we have A^min(l) = 3, so a minimum number of 3 gratings is required to form 
a grating echo. In this subsection three D = 1 TLI schemes will be discussed using the 



displacement diagrams 



m 



Figure [6] corresponds to a TLI configuration that is most frequently explored in previous 



works 



26|. In this configuration, standing wave is pulsed at time Ti = 0, T2 = T/2, before 



the revived density grating is probed at time T3 = T + r. Equations (fT^ and (HM require 
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FIG. 6: Diagram corresponding to an "equal time 3-pulse" TLI. 
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FIG. 7: Diagram corresponding to a "balanced 4-pulse" TLI. 
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FIG. 8: Diagram corresponding to an "elongated 4-pulse" TLI. 



Til = —1, 77-2 = 2. To be specific, we shall call this configuration an "equal time 3-pulse" 
sequence. 

Figure [7] corresponds to a "balanced 4-pulse" TLI configuration [9|], with the standing 
wave pulsed at Ti = , T2 = T/4 and T3 = 3T/4, before the revived density grat ing 
is probed at T4 = T + r. Similar timing sequences have been explored in Ref. 29|] in 
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three configurations, namely, the "fast echo", "slow echo" and the "stimulated echo". The 
"balanced 4-pulse" configuration is understood as a timing sequence such that the "fast 
echo" and the "stimulated echo" in ref. 29(] become degenerate. As a consequence, more 



than one paths contributes to the revived density grating. Equations f[T^ and f[T5]) require 
ni = 2n — 3, n2 = 4 — 3n, 713 = n , where n is an arbitrary integer. The contribution from 
different paths add up coherently in Eq. ([9]). 

Figure [8] corresponds to an "elongated 4-pulse" TLI [8j, with the standing wave pulsed 
at Ti = , T2 = and T3 = T — Tg, before the revived density grating is probed at 
T4 = T + r. The "elongation" corresponds to the timing parameters T » Ts. We shall 
follow the terminology in Ref. 29|] so that the "elongated 4-pulse" configuration is considered 
as an "elongated stimulated echo". Equations ffT^ and ffT^ require rii = —1, n2 = 1, 71,3 = 1. 

The calculation of the interferometry output for the three configurations in Figs. O [3, El 
is straightforward using the rules in Fig. [2l The amplitude and phase of the revived pq are 
summarized in Table [H where in the last row we also include a k- vector modification factor 
due to the perturbation from a quadratic potential along Gx. The latter will be explained 
in detail in the next section. 



F. A two-grating-vector example 

The examples in the last subsection have been restricted to one dimension. To complete 
the discussion on the application of the displacement diagram, in this subsection we discuss 
a grating echo Talbot-Lau interferometry sequence, that involves a manipulation of atomic 
coherence using standing wave pulses with the k-vectors span two dimensions. 

As a particular example, we consider a 2D standing wave setup, that involves two sets 
of counter-propagating traveling waves with the k-vectors orthogonal to each other [see 
Fig. El^b)]. In particular, we consider the traveling wave light fields labeled with -Ei, 2,3,4, 
with k-vectors specified as ki 2,3,4. 

According to Eq. (fT7|) . a minimum number of N=2+2=4 pulses is required to complete 
the TLI sequence in the two dimensional case. An example sequence is given in Fig. M^a) 
with a 2D displacement diagram. The displacement diagram is a 3D plot that involves 
two spatial dimensions (reciprocal phase space X-Y) and a time dimension (t). For clarity 
purpose, the displacement line after each vertex is drawn in a different color. The standing 
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FIG. 9: Diagram corresponding to a 2D TLI. (a) The (X,Y)-t displacement diagram, (b) The 
arrangement of the four traveling wave k-vectors. (c) The projection of the displacement diagram 
into the X-Y plane (with arrows to indicate the direction of propagation) . The four graphs on the 
right bottom gives the relevant standing wave k-vectors at t = 0, T/4, T/2, and T. For convenience 
of discussion in the main text, the four diagrams are labeled as (d), (e), (f), (g). 

wave k-vectors are colored according to the color of the displacement line after the scattering 



vertex. The projection of the displacement lines into the X-Y plane is shown in Fig. [9]^c) |40 |. 
The standing wave pulse sequence is given in Fig. [9](a) as well as in Fig. [9] d,e,f,g. In this 
particular example, the atoms are diffracted by two standing waves at time t = 0, with 
k-vectors Qi = k2 — ki = 2key and Q2 = ks — ki = k{ey^ + Gy). The pulse area are 61 and 
9[ respectively. The second pulse at t = T/4 is with a single grating vector Qi and a pulse 
area 62- The third pulse at t = T/2 is with a single grating vector Q2 = 2A;ex and a pulse 
area ^3. The grating echo induced at time t = T + r is due to a probe light from Ei mode 
that is Bragg-scattered to i?4 mode. Notice that according to Eq. f|T5l) we generally have 
xi = Qatt. By applying Eq. ([9]) to Fig. [9l it is straightforward to write down the expected 
TLI output at time T + r: 

E^{T + t) = igJi{2ei sm2uJkT)Ji{2e[ sin2cjfcr) 
J_2(2^2 sin 2uJk{3T/8 + r)) J_i(2e3 sin2cUfc(T/2 + r))e-'(^)', 
with ujk = ^ and Vk = —■ 
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III. TLE IN A QUADRATIC POTENTIAL 



In the last section we have introduced the diagrammatic technique to calculate a general 
TLI grating echo that involves free space matter-waves interacting with N-1 standing wave 
pulses and a probe pulse. In this section we extend the technique and consider matter-waves 
confined (or anti-confined) in a quadratic potential. 

A. The displacement diagram with a quadratic potential 

We shall assume that the perturbation from the quadratic potential is ignorable during 
the standing wave pulse durations. Thus we only need to modify the (A3) part in the model 
in Section Hi A[ 

^W) = ( V + • + E cos(Q. ■ f + (24) 
Z Z • 

For a general discussion, in Eq. u; is a rank two tensor and can be purely imaginary 
for an "anti-trap" . 

1. Basic displacement diagrams 

The propagation of Weyl functions in a harmonic potential is derived in Appendix A. We 
have the modified propagation rule [compare with Eq. ([71)], 

W{q, X, t) = W{q ■ cos cut + x ■ u; ■ sin cut, x ■ cos cut — q ■ cj""^ ■ sin tut, 0). (25) 

The modified propagation rule is summarized with the diagram in FiglTOlfa). together 
with the scattering diagram and probing diagram in FigJTOlfb) and (c). Figure [TO] gives the 
basic displacement diagrams for TLI in a harmonic trap (compare with Figl2]). 

2. The displacement diagram for an N-pulse TLI 

In Fig. [TT] we sketch the diagram for a general N-pulse TLI weakly perturbed by a 
harmonic potential. The TLI output can be expressed formally identical to Eq. ([9]) except 
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\(t) = \ coso)t-q o ^ sin to? 



FIG. 10: Modified propagation rules due to Eq. ([25l) 



X 




FIG. 11: Diagram giving (x(t), q(t)) in a particular path in a general N-pulse TLI under harmonic 
perturbation (blue curve). For comparison, the red curve gives a displacement diagram without 
harmonic perturbation. 

here qi and Xi need to be evaluated according to Fig. [TOT a). we have: 

E^(T^) = -z^7E tn ^n,(2^.sin^^)]iy(qi,ii,Tf). (26) 

Equation (126!) can be evaluated with the displacement diagram in Fig. [11] and the rules 
given by Fig. [TOl and is the second major result of this work. In what follows we discuss 
TLE due to a harmonic perturbation in a weak (tuT << 1), pulsed {{uj(t) ~ S(t)), and 
strong {ujT >> 1) regime according to Eq. fl2^ . 
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B. TLE due to weak harmonic perturbations 



1. The dephasing factor due to modified input displacement line 

We apply Eq. fl2Bl) to the "uniform and broad" atomic sample specified by Eq. f llOl) . The 
Weyl function in Eq. fl26|) is given by 

iy(qi,ii,rf ) = /(?(r-/3p)/(p)e^'ii----Prf3rf/3p ^^^^ 
= G(qi)F(ii-/5qi), 

where G{q) and F{x) are the Fourier transform of g{r) and /(p) respectively. 

The slope of the input displacement line in Fig. [11] is modified to be qi(Li;,T) 7^ 0. The 
grating echo given by Eq. fl26|) has a contrast reduction factor C{T) = G'(qi) according 
to Eq. (1271) . For example, for a Gaussian or a Lorentzian distribution of atomic density 
distribution of a length 21, the contrast reduction factor is given by: 

C {T^Qaussian 6 (Si ) / ^ /qq\ 
'-'V j Lorenzian ^ 

The reduction of the grating echo amplitude as given by Eq. fl28l) is due to dephasing, 
as it stems from the inhomogeneous phase broadening for matter-wave at different location 
in the harmonic potential. Here qi is evaluated up to second order in ujT for the three ID 
examples discussed in Section III El (using Fig. [10]) and are listed in the last row of Table [B 
An interesting observation is that qi has the same T dependence as those in the fourth row 
that gives the acceleration induced phase shift. Of particular interest is the n = 2 path in 
Fig. [71 From Table [T] we see gi = and the harmonic potential induced dephasing is greatly 
suppressed. In this case higher order terms need to be evaluated, and it's easy to find to 
the fourth order, qi = ^(cjT)^. 

The dephasing described by Eq. ( [28l) is in principle harmless if we know exactly where 
and how strong the harmonic perturbation is, so that we can slightly modify the probing 
k-vector Qtv to compensate for any "magnification" or "demagnification" of the k-vectors. 
For a more general potential that is anharmonic, the dephasing results in a reduction of 
interference fringe contrast that cannot be easily compensated. To evaluate the dephasing 
factor due to a weak and smooth potential, one may approximate the anharmonic potential 
with several local harmonic potentials. 
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FIG. 12: An example of shifted echo time due to a weak harmonic confinement. The blue line 
corresponds to the displacement diagram in free space, while the red line corresponds to the 
displacement diagram due to a harmonic potential. 

2. Echo profile shift 



Equation fl271) also predicts a P dependent shift of echo time. Connected with a typical 
experimental situation where a localized thermal atomic sample is expanded during a waiting 
time T^ait, we replace the notation P — > T^ait- If « 1 and T^^jau = then xi ~ Xi. For 
a large T^ait, the factor F{QiyT — qiTwait) in Eq- (I2ZD shifts the envelope of Eq. (12^ as a 
function of r by an amount (we assume Qn and qi along the same direction for simplicity.): 

'echo 



5T 



Q 



-T 



wait- 



N 



(29) 



As an example, the echo time shi: 



displacement diagram in Fig. [12] 44 1. 



t in the "equal-time 3-pulse" TLI is illustrated with the 



3. Recoil phase shift 

To finish the discussion of TLE in a weak harmonic potential, we consider the factor 

N-l 

" n Jn i'^di sin ii^)" in Eq. ([9]). As discussed in the last section, for a time setting {Ti = 

i=i 

0,T2,...,T} = T{si = 0,S2,.-.,siy = 1} that satisfies the echo condition Eq. f|T5l) . E'^ is a 
periodic function of total interrogation time T when there are no harmonic perturbations. 

With harmonic perturbation, "sin in Eq. (Q is replaced by "sin in Eq. fl2BI) . 

that in general breaks the linear relation between the recoil phases in Eq. ( fTSl) and T. 
However, if the total interrogation time T is small so that uT << 1, the linear relation is 
still retained Parallel to Eq. (I29l) . we shall consider the atomic sample released from a 
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FIG. 13: A schematic optical layout representing the setup in Lau's original discovery [ij]. L 
is an imaging lens at downstream of two gratings Gi and G2 that are subject to an incoherent 
illumination. The screen S is put at the focal plane of the lens where a periodic image is formed, 
with its contrast depending on the separation a between Gi and G2. 

local source with a very long T^^ait, so that Xj = + qiT^ait- The amount of recoil phase 
shift is given by: 

= ^VQ, ■ q^T^au, (30) 

with vq^ = the recoil velocity. This shift is systematic in recoil frequency measurements, 
particularly if the measurement is carried out in a trap 39|] . The sensitive dependence of the 
echo shift in Eq. fl29l) and the recoil phase shift in Eq. fl30|) may also be explored to measure 
small curvature of a potential including gravity gradients. Comparing with the method of 
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42j, monitoring the echo shift or 



using two interferometers separated by a large distance 
recoil phase shift of a time domain grating echo interferometer is more convenient if a field 
curvature over short distances is of concern. 



C. Pulsed harmonic trap and a time-domain Lau effect 

In this subsection we shall consider the modified propagation rules given by Eq. (!25l) in 
the short-pulse or the thin-lens regime. We then apply the result to study the Fraunhofer 
diffraction of incoherent illumination over two gratings (the original discovery of Lau [l4|). 

We consider a pulse of a one dimensional harmonic trapping potential with the trap 
frequency uj switched on for an interaction time r^, short enough so that uJTh << 1 and the 
atoms do not travel significantly in the harmonic trap during the interval (e.g., the Raman- 
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FIG. 14: displacement diagram for a time-domain Lau effect 



Nath regime). We may effectively consider uj(ty = fS(t), with / = oo^Th- Equation ( 125|) is 
reduced to 



The pulsed harmonic potential is considered as a thin lens in the time domain. To see 
this, consider a paraxial atomic beam with beam velocity Vbeam- If we make the replacement 
z := Vbeamt and fz ■= "''y'" , the effect of the pulsed harmonic trap described by Eq. (I3T!) to 
the atomic beam is equivalent to a lens with focal length 

As a particular application of Eq. fl31l) . in the following we consider a time-domain version 
of Lau's original experiment. We consider an optical layout sketched in Fig. [131 A un- 
collimated light passes through two gratings Gi and G2, with grating constant d, separated 
by a distance zi^2- The Fraunhofer diffraction is observed using a lens with focal length fz 
on a screen 'S'. As discussed in the introduction, a similar setup was used by Ernst Lau in 
1948 [isl to make the original observation of the Lau effect: Interference fringe for light 
at wavelength A appears on the screen 'S' when 2:1 2 = 

A similar experiment can be implemented in the time domain using standing wave pulses 
and a pulsed harmonic potential. To simplify the discussion, we consider the standing wave 
with k-vector Q = Qe^^ and the pulsed trapping potential is also along the direction. The 
standing wave is pulsed at time Ti = and T2 = Tg. The harmonic trapping potential is 
pulsed at time T at an amplitude u with a short duration Th and we have / = uj'^Th- The 
correspondent displacement diagram is shown in Fig. [TH We again consider a "uniform and 
broad" atomic sample [Eq. (fTOj) ] with coherence length Ac. From the diagram we have the 



W{q, X, 0+) = W{q - fx, X, 0"). 



(31) 
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revival of an n order grating: 

PnQ/M{t = T + 1// + r) = W{q = nQ/M, x = 0,t = T + 1/ f + r) = 
U2e, sin(^)) J_„(2^^2 smin'ujQiT, + r/M)))e-^(— ) , 

where M = -rir is a magnification factor that demagnifies the fringe k-vector from nQ to 
nQ/M so that the fringe spacing is magnified by M. 

We notice the similarity between the "elongated 4-pulse" configuration in Fig. [8] and 
the time-domain Lau configuration in Fig. [HI The time-domain Lau effect may be most 
easily realized with a pulsed optical dipole trap. By reducing Tg and the trap pulse strength 
/ = uj'^Th, the magnification factor M can be fairly large. In this way, wavelength-scaled 
fringes can be converted to a much larger size, resolvable using regular optical imaging 
system. 



D. TLE in a tight harmonic trap 

In this subsection we consider the revival of atomic density grating in a harmonic trap. 
Motion in a harmonic trap is periodic and one expects the grating echo to revive due to 
harmonic confinements. In particular, after a standing wave pulse the atoms bunch toward 
a standing wave potential minima to make a density grating, that repeats itself as the 
displacement line in Fig. [15] crosses the x=0 axis. Consider the initial atomic distribution 
to be a Maxwell distribution in the harmonic trap with spatial extension L and coherence 
length Ac. We see that in Fig. [T5]r2i = ±1, and expect a TLI output: 

E\T) = -«(?Ji(2^isin(^^sinu;r))(e-«i-'=°^-^)«^)' - e-((i+™«-^)Q^)')e-(^^'^("^»', (33) 

u 

with VQ = ^,and UQ = ^. 

In Fig. [16] we plot examples of the atomic grating revivals. The trap frequency uj is chosen 
to be 0.14 times the recoil frequency uoq, while the initial atomic phase-space distribution 
is tuned via temperature. Figure [TBT a) corresponds to the atoms initially occupying the 
ground state of the harmonic trap only. Figure [TBT b) corresponds to the initial condition of 
a thermal distribution with temperature equal to 50 times the recoil temperature. In both 
cases the trap frequency determines the period of the grating revival. 

Complimentary to the proposal of using the echo shift [Eq. (1291) ] or recoil phase shift 
[Eq. ([30]) ] to measure very weak harmonic perturbations in the last section, the revival 
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FIG. 15: Diagram of the single pulse TLX in a harmonic trap 




FIG. 16: Representative revived echo signal calculated according to Eq. ([33j) . The trap frequency 
CO is chosen to be 0.14 times the recoil frequency, 9i = 2. Figure (a) corresponds to an initial 
atomic sample at the ground state of the harmonic potential. Figure (b) corresponds to the initial 
condition of a thermal distribution with temperature equal to 50 times the recoil temperature. 

of the atomic density grating predicted by Eq. (!33l) may be experimentally exploited to 
precisely determine the trapping frequency of tight traps. 

IV. DISCUSSION AND SUMMARY 

A. Grating with general transmission functions 

The discussion so far has been restricted to grating potentials using sinusoidal light shifts. 
In the following we shall discuss the generalization of the formula to include arbitrary grating 
transmission functions. 

We shall first modify the light shift in Eq. ([3]) to include an imaginary part, e.g., Qi+iTi = 
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6i5{t — Ti), with 9i = 9f + iO^. The imaginary part of the hght shift potential could be due 



to the optical pumping effect of the light pulses 



461] that act as an amplitude mask for 



could also describe a particular harmonic 



16j. Following the principle of analytical 



atoms. More generally, the imaginary potentia_ 
component of an absorptive material grating 
continuation, Equation ([H]) and the scattering diagram in Fig. Ulp) can be modified to 
include this imaginary part as (Appendix A): 

W{q, X, Tt) = E nQ^, x)iy(q - nQ„ x, Tr), (34) 

n 

where 

Fie,, nQ„ x) = e-'''^,UV^^^), 
a = 2ef sin b = 291 cos 
For a general amplitude/phase grating that has more than one spatial Fourier component 
in its transmission function, in Eq. (IMI) F(6'j, nQi, x) has to sum over contribution from all 
the Fourier components for the i*^ grating, e.g.. 



F(^„nQ„x)= e-'''^\r-^JU^M^f), (36) 

Is. \ CLl + bi ^ 

ly.m=n v i i i 

where 6i{l) give the effective pulse area due to the Z*^ order harmonic of the grating potential, 
ai = 26^(1) sin ^ and 6/ = 2^/ (/) cos 

Notice here in Eqs. ( IMl) . ( l35i) . ( l36l) and also in Eq. (IHD the subindex i is used to refer to 
the i^'^ pulse in a TLI in Fig. [3], and can in fact be ignored for convenient reading. 

To detect the resulting atomic density grating as the output of TLI, in this paper we 
have considered the method of a grating echo using Bragg scattering of a probe light, that 
is formulated in Eq. ([5]). Another detection method frequently used for matter gratings and 
for resonant standing wave masks is instead to monitor the overall transmission after an 



amplitude grating [l6|, ll8|, l25|, l30|, l3l|]- The formula developed in this work, e.g., Eq. ([5]) 
and Fig. [2](c) can be modified accordingly to (coherently) sum over the interaction with all 
the spatial harmonics of the transmission grating. We end up with the transmission signal 
S(Tn) given by [compare with Eq. ([5])]: 

S{Tr,) oc E A; / dhp{r, T^)e^lQ-•^ (37) 
I ■' 

where Ai is the coefficient of the order spatial harmonic of the N^'^ grating transmission 
function. 
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B. Matter- wave coherence, displacement diagram and recoil diagram 



In this work we discussed using Weyl functions to study the Talbot-Lau effect (in three 
dimensions) in free space and in a quadratic potential. We now come back to Eq. ([T]) 
to discuss the Weyl function as a measure of matter-wave coherence, and compare the 
"displacement diagram" technique here with the "recoil diagram" introduced by the authors 



of ref. 
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that has been widely used in the atom interferometry community [1 



1. Displacement diagram and conservation of matter-wave coherence 

We come back to Eq. ([T]). As mentioned in the introduction, lV(q, x, t) measures the over- 
lap, or the second order field-correlation of matter- waves displaced in phase space by (/iq, x). 
As a particular example, W{(\ = 0,x, t) measure the spatial coherence of a wavepacket be- 
fore and after a particular displacement x, which is also called the longitudinal coherence 



in matter- wave interferometry 48|. In 48|] it was pointed out that longitudinal coher- 



ence is conserved during free evolution, which is a fairly counterintuitive observation since 
wavepackets of atoms may expand a lot and one may suspect the coherence length of the 
wavepackets increases accordingly. 

The conservation of W{c{ = 0,x, t) during free evolution is only a particular example 
of Eq. ([7]), which more generally describes the conservation of matter- wave coherence that 
propagates in reciprocal phase space (q, x). In particular, momentum coherence of matter- 
waves given by W{c{, x = 0) = Pg is a Fourier component of atomic density grating (fringes), 
that propagates according to iy(q, x = — qt) but with conserved magnitude. These propa- 
gating coherences reflect the correlation of matter-waves at different locations with different 
velocities, though they cannot be measured directly with any local measurement, except 
those on the axis x = 0. Equation ([8]) or flM|) show how these phase space coherence can be 
manipulated with periodic potentials to create a revived momentum coherence W^(q, x = 0) 
for a local measurement. 

Even more generally, it is easy to see from Eq. ([T]) the identity 

y"|H^(q,x,t)|2rfqrfx = Tr(p^), (38) 



which is equal to < p{t) >, the purity or linear entropy of matter- waves 49|,|50|. Matter- wave 



entropy is conserved in a non-dissipative/isolated potential, which may also be interpreted 
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FIG. 17: a) Recoil diagram of matter-waves due to a standing wave pulse. Circles on each line 
are "billiard balls" with radius equal to the coherence length of wavepackets (the correspondent 
displacement diagram is given by Fig. H]). b) Recoil diagram of a 3-pulse grating echo interferometer 
(the correspondent displacement diagram is given by Fig. [6]) . Standing wave is pulsed at time t = 0, 
T/2. At around time T the billiard balls along each line cross each other to create interference. 
The circles of "billiard balls" are not shown. 

as a conservation of matter-wave coherence through the Weyl function, given by: 

d r 

— / |iy(q,x,t)prfqrfx = 0. (39) 



2. Displacement diagram and recoil diagram 

The recoil diagram was introduced by the authors of ref. [3(| together with a "Bilhard 
Ball Model" of atom wavepackets to explain the photon echo phenomena, and also the 



Ramsey interference and two-photon echo in matter-wave interferometry [47|. In a (x-t) 
recoil diagram such as Fig. [T7l(a). lines represent the trajectories of centers of wavepackets 
diffracted by optical pulses. The wavepackets are represented with "billiard balls" with 
a fixed radius equal to the coherence length of the wavepackets. Matter-wave interferes 
whenever two billiard balls cross each other. The visibility of the interfering fringe relies on 
the overlap of the billiard balls, while the fringe position is determined by the relative phase 
between the two wavepackets. A rigorous justification of usin g th e "billiard ball model" 
is presented in which is related to the conclusion in ref. [48|] and also our previous 
discussion. 

Examples of recoil diagrams due to a standing wave pulse (displacement diagram in 
Fig. H]), and for an "equal time 3-pulse" (displacement diagram in Fig. [6]) grating echo 
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interferometer is shown in Fig. [T7] (a), (b) respectively. In contrast to the recoil diagram 
that concerns the coherence between individual pairs of wavepackets [such as the green path 
in Fig.fTWb)]. the displacement diagram describes the overall coherence of the atomic sample, 
which is a sum of the coherence amplitudes between all pairs of the wavepackets with a same 
relative phase-space displacement. The displacement diagram is particularly convenient to 
describe the overall coherence of multiple diffraction orders, particularly when a sequence 
of multiple diffractions is considered (such as an N-pulse sequence given by Fig. [3]). The 
overall coherence of an atomic sample is directly related to interference readouts in most 
interferometry setups where the interference fringes are retrieved across the whole atomic 
sample. 



3. Dephasing and decoherence 

Finally, it is instructive to consider the dephasing effect in both the recoil diagram picture 
and in the displacement diagram picture. For a matter-wave interferometer perturbed by 
an external potential, the phase of wavepackets may be in-homogeneously broadened that 
reduce the overall coherences of matter- waves and thus the fringe contrast. This is referred 
to as a dephasing effect. Dephasing effects can be described by a recoil diagram, by assigning 
each "billiard ball' with a different phase and by varying the "radius of the billiard ball" to 
account for wavepacket distortions. 

In the displacement diagram picture, the dephasing effect broadens and shifts the Weyl 
function distribution in the reciprocal phase space, for example, the dephasing effect due to a 
quadratic potential corresponds to a shift of Weyl function described by Eq. fl27j) . However, 
any broadenings or shifts of Weyl function distribution need to obey the conservation law 
given by Eq. fl39l) . In contrast, the reduction of interferometer fringe visibility accompanied 
by a violation of Eq. ( l39l) is due to a loss of information or an increase of entropy in the 



matter-wave system, and should be referred to as decoherence effects 
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52| 



C. Summary and conclusion 

We have generally discussed a vector theory of Talbot-Lau effects due to N standing 
wave pulses by considering an N-pulse grating echo interferometry configuration (Fig. [1]). A 
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displacement diagram (Figs. [2] and [TO]) is introduced to calculate the grating echo diagram- 
matically. To illustrate the convenience of the diagrammatic technique, we discuss examples 
of TLI that are related to recent experiments 8|, |9|, |39| , in particular, we discuss an example 
TLI with standing wave k-vector span two dimensions. We study the dephasing of these 
TLI schemes due to harmonic confinements (Table [I]), and predict the echo shift [Eq. (1291) ] 
and the recoil phase shift [Eq. (l30l) ] that have been observed experimentally [10]. These 
effects should be useful to measure the curvature of an atomic potential. The displacement 
diagram method is further extended to include a pulsed harmonic trap [Eq. (13T1) ] and we 
discuss a time-domain version of the original Lau effect [Eq. (l32l)]. Using the displacement 
diagram in a tightly confined trap, we give the expression of the atomic density grating 
revivals [Eq. fl5^ ]. The formula introduced in this work can be extended to describe TLE 
due to diffraction gratings with general grating transmission functions. Application of Weyl 
functions in this work has been restricted to Talbot-Lau effects involving multiple diffrac- 
tion orders from diffractive gratings. But the method is generalizable to other interferometry 
configurations with fewer interfering paths. 

Coherences, or "off-diagonal terms of a density matrix" ^], have a clear definition in 
discrete systems but not so clear for matter- waves that have a continuous spectrum js^]. In 
the last part of this work we have briefly discussed Weyl function as a convenient measure 
of matter-wave coherence. Weyl function as a phase-space correlation function measures 
the overlap of matter- waves before and after a phase-space displacement. Similar to the 
definition of the Wigner function. Equation ([T]) involves a Fourier transform of off-diagonal 
part of a density matrix. The Fourier transform removes the ambiguity due to the ill-defined 
off-diagonal terms, while keeping all the physical information of matter-wave coherences 53). 
We thus suggest using Weyl function to define matter-wave coherence and more generally 
to define coherences in system with a continuous spectrum. We suggest that according to 
this definition, a conservation law of matter-wave coherence [Eq. ( !39l) ] that directly link 
to the matte-wave entropy can unambiguously distinguishes decoherence and dephasing 
effects. Practically, it seems that to consider the Weyl function or the displacement operator, 
instead of considering the evolution of particular matter-wave states, is more convenient 
for elucidating the dynamics of matter-wave coherences during a particular interferometric 
observation. 
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APPENDIX A: PROPAGATION PROPERTIES OF THE DISPLACEMENT OP- 
ERATOR 

We consider the evohition of the displacement operator D{c{, x) — e'*^*^'"""^ ^^ in the Heisen- 
berg picture. Notice that D{q,x) can also be expressed as e''^^e^^'^e~^^ . 
We consider three type of propagators: 



~2 

3— — |-a-f)t 



ft2 ,,2-2, 



Uf{t) = e 

Uharmonicit) = e-(^+^)*, (A-1) 
TT „— iScosQ-f 



We have, during free propagation in an accelerating frame with acceleration constant a, 

(A-2) 



D(q,x,i) = Uj\t)e'^'i-^-''-P^Uf{t) = 



In a harmonic trap: 

giq-(f-cosu)t— p-u)~^-sinu;t)+ix-(p-cosajt+f-w-sina)t) _ (A-3) 

i)(q • cosa;t + x • a; • sina;t,x • cosa;^ — q • ■ sina;i). 
Transformation of ^)(q, x) under a standing wave pulse is given by: 
D+(q,x) = t/,;ie-^e^^-e-^C/,,Q = 

g-i^giqfgiecosQ-(f+f )g-t(9cosQ-(f-f ^ (A-4) 

g_i2^giq.fgi2esin(%i)sin(Q.f)g-«^^ = ^ (2^ sin (q + hQ, x) . 

n 

Finally, we consider the propagator due to a standing wave pulse with complex area: 

^' ' (A-5) 

^Q,e — 6 
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The transformation of £)(q, x) is given by: 



Ul^e 2 e 



XJJ 



e 2 
e 



,Q 



(A-6) 



E E ^n^+n(2^^ sin ^)/^(2^^ COS ^)^-(q + nQ, x). 



In Eq. flA-61) /„(x) = z "J„(ia;) is the modified Bessel function of the first kind. From 
Eq. (1A-6P to Eq. we use the Graf summation formula for Bessel functions with complex 



arguments 551]. The multiple- va^ 



be single valued according to [55|. 



ued arguments of Bessel functions in Eq. (l34l) is chosen to 
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